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The physical origin of the strong CP problem in QCD, rooted in the structures of the vacuum of 
the standard model, is reviewed. The chiral solution to this problem, with its accompanying axion, 
is explained and various characteristics of axions are detailed. Although visible axion models are 
excluded experimentally, models with very light and very weakly coupled axions are still tenable. 
The astrophysical and cosmological implications of such axion models are discussed, along with 
ongoing experimental attempts to detect such, so called, invisible axions. 
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Q\ • THE U(1) A PROBLEM 
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Quantum Chromodynamics (QCD) describes the strong interactions of hadrons in terms of the interactions of their 
C , quark and gluon constituents. The QCD Lagrangian 
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7— I ' for / flavors of quarks has a large global symmetry in the limit when irif — > : G = J7(/)r x E/(/)l- This 
symmetry corresponds to the freedom of arbitrary chiral rotations of the / flavor of quarks into each other. Because 
1 m u ,md <C Aqcd — with Aqcd being the dynamical scale of the theory — in practice only a chiral U(2)n x f (2)l 
symmetry is actually a very good approximate global symmetry of the strong interactions. This symmetry, however, 
is not manifest in the spectrum of hadrons. It is instructive to understand why this is so. 

The C/(l)v subgroup of the £/(2) R x U(2) L symmetry, corresponding to vectorial (V = R+L) baryon number, is an 
\Q , exact symmetry of QCD, irrespective of the value of the quark masses. In the limit of m u — rrid a further SU(2)v 
0^ ' subgroup would be an exact (isospin) symmetry of QCD. Because the u- and d-quark masses are so light compared 
to Aqcd, even with m u ^ one has an approximately degenerate multiplet of hadrons corresponding to this 
symmetry (e.g. the 7r-triplet and the doublet of neutrons and protons). However, the corresponding axial symmetries 
in J7(2)r x ?7(2)l (A = R — L) are not seen in the spectrum. For instance, there are no parity doublets degenerate 
with the neutron and proton. This phenomena is understood because the SU(2)a symmetry is not preserved by the 
QCD vacuum. Indeed, as is appropriate for a spontaneously broken global symmetry, there is instead in the QCD 
hadron spectrum an approximate triplet of Nambu-Goldstone bosons — the pions — whose mass vanishes in the limit 
that m u ,m<i — > 0. Remarkably, however, there is no corresponding pseudoscalar state with vanishing mass, in the 
same limit, corresponding to the Nambu-Goldstone boson of a U{1)a symmetry. 

The nature of the approximate SU(2\y x 5f (2)a x f/(l)v symmetry of the strong interactions was actually 
understood before the advent of QCD m. Although the anomalous features of the f (1)a symmetry were also 
pinpointed around the time of the development of QCD as the theory of the strong interactions |2|, the resolution 
of the U(1)a problem required a better understanding of the QCD vacuum. In short, the reason why there are no 
approximate Nambu-Goldstone bosons associated with this Abelian chiral symmetry is that, as a result of chiral 
anomaly this is really not a quantum symmetry of QCD. In turn, however, this more complete understanding 
opens up another problem — the strong CP problem Namely, the richer vacuum of QCD, combined with the 
violation of CP in the weak interactions, allows for the presence of an effective interaction 

— C^s ~ 

^strong CP = GauvG^ (2) 
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which leads to an enormous neutron electric dipole moment \d n ~ ed(m q /M^)^ , unless the parameter 9 is tiny 
((9<I(r 9 ). 

In QCD, the dynamical formation of quark condensates, (uu) — (dd) ^ 0, breaks the U(2)a global symmetry. As a 
result, in the limit of vanishing quark masses, the pion triplet are the Nambu Goldstone bosons associated with the 
spontaneously broken SU(2)a current = 97^75 Qq, where q is the (it, d) doublet of quark fields. The existence of 
these massless states produces a q 2 = pole in the Green's functions of the SU(2)a currents with pseudoscalar quark 
densities. One has 
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Mtj(q) = J d 4 xe-^(T(j£(x),q l5 r iq )) 



(qq)5 tJ q' 1 

~ 5= h non singular terms . (3) 
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Because the (qq) condensates also break the ?7(1)a global symmetry, naively one expects also a similar q 2 singularity 
associated with the U(1)a current, connected to a further Nambu-Goldstone boson-the r\ in this 2-flavor discussion. 
However, such an approximate Nambu-Goldstone excitation cannot really exist, since if it did, it would lead (in the 
limit of finite but small quark masses) to a state essentially degenerate with the pions Q — in contradiction with 
experiment. 

The situation is a bit more complicated because the U(1)a current J§ = 97^755? has a chiral anomaly [|| 

« = ^G^G a , v . (4) 
However, since the RHS of the above equation is a total divergence |5j 

(5) 



one can construct a conserved chiral current 



Gaf3-y T^fabc ^-cf 
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As a result, in the massless quark limit, if (qq) 7^ 0, one can show that the Green's functions of the j£ current have 
a. q 2 = singularity analogous to that of the chiral SU(2)a currents: 



M»(q) = J d 4 xe~^(T(4(x), q l5 q)) 
(qq)q^ 1 



1 non-singular terms . (7) 



The C/(1)a problem, then, is why such a singularity is not manifested physically in the spectrum of hadrons? 



THE QCD VACUUM 

The resolution of the C^(1)a problem came from a better understanding of the structure of the QCD vacuum. For 
gauge fields one can take the vacuum state to be the state where the vector potential is either zero or is in a gauge 
configuration equivalent to zero. A proper study of these configurations reveals that the structure of the vacuum 
state for non-Abelian gauge theories (like QCD) is richer than expected. Following Callan, Dashen, and Gross Q it 
is convenient to study these theories in the temporal gauge A° — and concentrate on an SU (2) subgroup. With this 
gauge choice, the spatial gauge fields are time independent and under a gauge transformation transform as 

^AUr) = A l (r) -> nmAH^VLif)- 1 + —0(0 Vfi -1 ^ ■ (8) 
2 g s 

Thus, in this gauge, vacuum configurations corresponds to a vanishing vector potential, or to — fi(r)V 4 f2 (f). 

The rich vacuum structure of QCD results from the requirement that the gauge transformation matrices £l(r) go 
to unity at spatial infinity. Such a requirement maps the physical space onto the group space and this S3 — > S3 
map splits O(r) into different homotopy classes {f2 n (r)}, characterized by an integer winding number n detailing how 
precisely f2(r) behaves as r — > 00: 

il n (r) ^ e 2m . (9) 

r — * 00 

The winding number n is related to the Jacobian of the S3 — > S3 transformation 

n= i&J ^ ^ A n(Mi{r)A k n {r) , (10) 
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with A l n (f) being the pure gauge field corresponding to the gauge transformation matrix Cl n (f). 

Because one can construct the n-gauge transformation matrix r2„(r) by compounding n times f2i(r), it follows that 
the vacuum state corresponding to A n (r) is not really gauge invariant. Indeed, the action of the gauge transformation 
matrix f^i on an n-vacuum state gives an (n + l)-vacuum state 



fill In) = \n + l) 



(11) 



Nevertheless, one can construct a gauge-invariant vacuum state — the, so called, 0-vacuum-by superposing these n- 
vacua [pi 
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Indeed, it is easy to check that £li\0) = e \6). 

Because of this more complex vacuum structure, the vacuum functional for the theory splits into distinct sectors. 



+ {0\0)- = ^2^ m8 e- ine + { 
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This superposition of transition amplitudes from t = — oo to t = +oo with fixed n-vacuum difference v, weighted by 
different phases e w6 introduces an additional intrinsic source of CP violation. By using Eqs. (5) and (10), one can 
show |i| that the integer v has a gauge invariant meaning 



Thus, writing for the vacuum functional the usual path integral formula 



(15) 



(16) 



one can re-interpret the 6 term in terms of an effective action 

S cS = S[A] + 6^- I d A xG^G atlv 



(17) 



This additional term violates P and CP, since it corresponds to an E a ■ B a interaction of the color fields. 

Perturbation theory is connected to the v = sectors where the GG term vanishes. The effects of the v ^ sectors 
thus necessarily are nonperturbative. These contributions are naturally selected by the connection of GG to the chiral 
anomaly. For nt flavors, the axial current j£ has an anomaly 



8 Jt* — nt-ln fiw 



Thus a chirality change AQ5 is related to v as: 



AQ 5 = / d"xd^ = n f -^- I d^xG^Ga^ =n f u 



(18) 



(19) 



The solution to the U(1)a problem is connected to the chirality breakdown in QCD in the v 7^ sectors. If these 
sectors are included — as they should — then U(1)a is never a symmetry. 

For rif = 2, the conserved current discussed earlier is not gauge invariant for transformations with non-trivial 
winding number. One can show Q that the associated chiral charge 



d 3 xJ° 



(20) 



transforms as 
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fiiQs^r 1 = Qs + 2 (21) 

when one worries about v ^ sectors, 't Hooft showed that it is this lack of gauge invariance which is crucial to resolve 
the Z7(1)a problem 0. Although Green's functions of the J§ currents contain 1/q 2 singularities, these singularities 
are cancelled for the corresponding Green's functions of the gauge invariant currents j£. The Green's function for 
this current with a pseudoscalar density obeys the equation 

M"(«) = M"(q) + d 4 xe-^(T(K»(x), q l5 q)) (22) 

Although the first term above is singular at q 2 = 0, when one considers the effects of v + sectors, so is the second 
term. These singularities cancel for the Green's function of the physical current J§ |7J. Because of Eq. (19), the 
C/(1)a symmetry is really not a symmetry at all. Therefore, not surprisingly, there are really no physical Nambu 
Goldstone bosons associated with this pseudo symmetry! 

THE STRONG CP PROBLEM 

The effective interaction (17) is not the only new source of CP violation arising from the more complex structure of 
the QCD vacuum. It is augmented by an analogous term coming from the electroweak sector of the theory. The origin 
of this additional interaction is as follows. In general, the mass matrix of quarks which emerges from the spontaneous 
breakdown of the electroweak gauge symmetry is neither Hermitian nor diagonal: 

Anass = </H U,„</|. " qu{M^)ijq Kj . (23) 

This matrix can be diagonalized by separate unitary transformations of the chiral quark fields. These transformations 
encompass a chiral J7(1)a transformation 

«a -» e M/2 g R ; <zl -» e~ W V (24) 
with a = — Arg det M. Such chiral transformations, in effect, change the vacuum ang le 8 . That this is so follows 

nf 1 — v 

from Eq. (21). Using this equation (replacing the factor of 2 by the number of flavors n/), it is easy to see that a 
chiral E/(1)a rotation on the 9 vacuum shifts the vacuum angle by Arg det M B: 



iaQ 5 



= \6 + n f a) = \6 + Arg det M) . (25) 



Hence, in the full theory, the effective CP-violating interaction ensuing from the more complex structure of the QCD 
vacuum is that given in Eq. (2), with 

9 = 9 + Arg det M . (26) 

The strong CP problem is really why the combination of QCD and electroweak parameters which make up 9 should 
be so small. As we remarked upon earlier, the effective interaction (2) gives a direct contribution to the electric dipolc 
moment of the neutron and the strong experimental bound on d n requires that 9 < 10 -9 In principle, because 9 
is a free parameter in the theory, any value of 9 is equally likely. However, one would like to obtain an understanding 
from some underlying physics of why this number is so small. 

THE CHIRAL SOLUTION TO THE STRONG CP PROBLEM 

Almost twenty years ago, Helen Quinn and I suggested a dynamical solution to the strong CP problem. What 
we postulated was that the full Lagrangian of the standard model was invariant under an additional global chiral 
U(l) symmetry. If this J7(1)pq symmetry existed, and it were exact, the strong CP problem would be trivially solved 
since 9 could be set to zero through such a chiral transformation (cf. Eq. (25)). However, physically such a symmetry 
cannot be exact. Nevertheless, what Quinn and I jl0| showed was that, even if J7(l)pq is spontaneously broken, the 
parameter 9 is dynamically driven to zero. However, in this case there is an associated pseudo Nambu- Goldstone 
boson in the theory, the axion |ll|].F] 



1 The axion is not massless because the chiral J7(1)pq symmetry is anomalous. As a result, the axion gets a mass of order 
m a ~ Aq CD //, where / is the scale associated with the breakdown of the C/(1)pq symmetry. 
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It is useful to understand schematically how a global Z7(1)pq symmetry solves the strong CP problem jj]. Basically, 
what happens is that by incorporating this symmetry in the theory one replaces the static CP violating parameter 
9 by the dynamical (CP conserving) interactions of the axion field a(x). Because the axion field is the Nambu- 
Goldstone boson associated with the spontaneously broken £/(1)pq symmetry, this field translates under a C/(1)pq 
transformation. If a is the phase parameter of this transformation and / is the scale associated with the breakdown 
of the symmetry, one has 

a(x) -> a(x)+af. (27) 

It follows, therefore, that if the effective Lagrangian describing the full theory is to be £/(1)pq invariant, the axion 
field must only enter derivatively coupled. This statement would be exactly true if the Z7(1)pq symmetry were not 
anomalous. Because of the chiral anomaly, however, C e s must also have a term in which the axion field couples 
directly to the gluon density GG, so as to guarantee that Jpg has the correct QCD chiral anomaly. 

The above considerations fix the form of the effective standard model Lagrangian, if it is augmented by an extra 



C/(1)pq global symmetry as we suggested 10 . Focusing only on the extra terms involving the axion field one has 



^sm — Am + Q-^—G% u G a uv — -d^ad^a 
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Ant. [d»a/f; V>] + jd^G^Ga^ . (28) 



In the above ip stands for any field in the theory and £ is a model-dependent parameter associated with the chiral 
anomaly of the C/(1)pq current 



cv 

«q = ^GTGV • (29) 

The presence of the last term in (28) provides an effective potential for the axion field. Thus its vacuum expectation 
is no longer arbitrary. Indeed the minimum of this potential determines the axion field VEV (a) . One has 



\ da I f8n [ a afn " 



= . (30) 

(a) 



What Quinn and I showed |I(J is that the periodicity of the pseudoscalar density expectation value (GG) in the 
relevant ^-parameter, 9 + forces the axion VEV to take the value 

(a) = - f j9 . (31) 

This solves the strong CP problem, since £§M' when expressed in terms of the physical axion field, a p h ys = a — (a), 
no longer contains the CP violating 9GG term. Furthermore, expanding V c g at its minimum, one sees that the axion 
itself gets a mass: 



a _ / d 2 V cS \ _ tja s d - , 

a ~\ 9a 2 / f87rda [ a 



(32) 

(a) 



Therefore, the standard model with an additional C/(1)pq symmetry no longer has a dangerous CP violating interac- 
tion. Instead, it contains additional interactions of a massive axion field both with matter and gluons characterized 
by a scale /: 

eff 1 

Asm = Am + Ant.^aphys.//; V 1 ] — 2^ a P h y s -^ a phys. 



2 m a a phys. + Pi J S ' t-^G% V G alJiV . (33) 
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AXION DYNAMICS 



In practical applications, it is more convenient to use the freedom of [/(1)pq transformations to replace the last 
term in Eq. (33) in favor of effective interactions of axions with the light pseudoscalar mesons (it and 77) of QCD 
p2| . Although the properties of axions are model-dependent, this dependence can be isolated in a few numerical 
coefficients related to the axion mass, its coupling to two photons and the mixing of axions with the neutral pion 
and the eta. In terms of the pion decay constant ~ 92 MeV and the SU(2) x U(l) VEV v ~ 250 GeV, it proves 
convenient to define a "standard" axion mass parameter 

mg = m ' f \ ^% * 25 KeV. (34) 
v (m u + ma) 

Then, for all models, one can characterize the mass of the axion and its it and rj couplings by 



and 



m a = X m mf ( - ) (35) 



with A m , A3, and A Q model parameters of O(l)J^) Similarly, one can characterize the effective couplings of axions to 
two photons by the interaction 



r - —K ° phys - v^F (V7\ 

47T / 



where K ail is again a model-dependent parameter of 0(1). It is clear that if f n <C /, the axion is both very light 
and very weakly coupled! 

It is informative to sketch the derivation of the above formulas in the original Peccei-Quinn model jl0|. Here the 
U(1)pq symmetry is introduced in the standard model by having two Higgs doublets: 

^Yukawa = I^Gl^IUR, + I ^l. <M| { + h.C (38) 

The presence of $1 and $2 allows £yukawa to be invariant under independent rotations of ur and c?R-the desired 
chiral U(1)pq transformation. The axion in this model is the common phase field of <!>! and $2 orthogonal to weak 
hypercharge. If x — v^jv-y is the ratio of the Higgs VEV and / — \/ v\ + v\ = {\f2Gp) — v ~ 250 GeV, then it is 
easy to isolate the axion content in <I>i and $2 as 
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ixa/f . <jj J^2_ 

' 2 V2 



e la/xf . (39) 



Using Eq. (27), to guarantee invariance of /^Yukawa, ur and d R transform under a J7(1)pq transformation as 

—tax 1 — ia/xi f An\ 

u R —> e mr ; d R -> e 1 d R . (40) 
Whence, it is easy to see that the current associated with the C/(1)pq symmetry isj^| 



y 1 - ■ y 



X 

i=l i=l 

i— 1 



2 In certain models, it is possible for A3 or A to vanish identically. 

3 Here I assumed that the charged leptons £r have a Yukawa coupling involving $2 . 
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This J7(1)pq current has both a QCD and an electromagnetic anomaly 



-OL s 



where [Q 



£7=^3 (*+- 



(42) 



(43) 



The QCD anomaly renders a bit more difficult the calculation of the effective interactions of axions with light 
hadrons [ fl3| |. Perhaps the simplest way to proceed is by using an effective Lagrangian technique to describe the 
interactions of pions and the r), both among themselves and with axions. Since ir and r] are the Nambu-Goldstone 
bosons associated with the approximate global C/(2)r x £/(2)l symmetry of QCD, their interactions are described by 
an effective chiral Lagrangian 



C 



chiral 



-/ 2 Tr c^E^E 



(44) 



with 



E = exp — [r • 7? + rj\ 

Jit 



(45) 



This Lagrangian must be augmented by a mass breaking term that mimics the Yukawa interactions (38) and thus 
involves the axion field, along with an axion kinetic energy term 



where 



A = 



-ixa I f 







D -ia/xf 



M 



m u +m d 









(46) 



(47) 



The presence of A guarantees the U (1)pq invariance of £ ax ion since under U (1)pq, to mimic the quark transformations, 









ia/x 



(48) 



The anomaly interactions, which break the above U(2) x U(2) x C/(1)pq symmetry through the coupling of gluons to 
axions and the rj, serve to give an effective mass term to the field combination which couples to GG. One hasQ 



-'anomaly 









■=)•] 







(49) 



with m 2 ~ m 2 S> m\ 



The quadratic terms in £axion along with £ a nomaiyj when diagonalized, allows one immediately to compute the 
axion mass in the model, as well as the axion mixings with tt° and 77 



O, — flphys. ^aTT^phys. ^arjVphys. 

A simple calculation yields for the Peccei-Quinn model p0[: 



(50) 



4 The factor of (N g — 1) in £ an omaiy enters since here one must include only the contribution of the heavy quarks in the 
axion-gluon interactions, as the light quark interactions are explicitly accounted for in £ ax ion. 
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x ) ' ' 9 \~ ' x J {m d + m u ) 



X o = ^(l-N g )[x+-) . (51) 



Using these results, one can deduce readily the effective axion to two-photon coupling. The electromagnetic anomaly 
of the ir° and -q fields 



C = — 



tt° + 5rj_ 



F^F^ , (52) 



through the mixing with axions, gives an effective 077 coupling coming from the light quark sector. To this coupling, 
one must add the heavy quark contributiion. This is given by Eq. (42), except that for the parameter £ 7 one should 
replace N g by N g — 1. Adding these contributions together yields finally, for the model, 

Ka-fj = N g (x + -) mu . (53) 
V x m u +m d 



THE DEMISE OF VISIBLE AXION MODELS 

By changing the detailed way in which the Higgs fields $1 and $2 couple to the quarks and leptons, one can obtain 
a variety of axion models ||l4| . All of these variant models, as well as the original Peccei-Quinn model, have been ruled 
out experimentally. I will not discuss in detail here the different evidence against these weak scale (/ = v) axions, as 
this is reviewed in depth in Ref. 0. I will, however, give a few examples to give a flavor both of the models and of 
the experiments used to rule these axions out. 

If / = v, then the resulting axions are quite light, since m a — 25A m KeV. Unless A m 1, then m a < 2m e and 
these axions are also very long lived, since they can only decay via the a — > 27 process. These light, long lived, weak 
scale axions are ruled out by the non-observation of the process K + — > ir + a. One predicts [O] 

BR(K + -> 7T+a) ~3x 10" 5 A^ (54) 

while, experimentally, from KEK there is a bound [ jl5| |: 

BR(K + -> 7T+ nothing) < 3.8 x 1(T 8 . (55) 

One can avoid this bound in models where A = Oj^] However, experiments looking for axion production via nuclear 
de-excitation in reactors are sensitive to both A and A3. The absence of a tell-tale a — > 27 signal downstream, 
although less stringent than the K + — > ir + a process, then serves to rule out this variant [ fl6| |. 

Weak scale axions can be short-lived if m a > 2m e . For this to happen, A m must be large and this necessitates either 
x or x~ l to be large. This possibility runs into difficulty with experiments looking for axions in quarkonium decays, 
QQ — > 07. If the C/(1)pq assignments of the Higgs fields are as in the original Peccei-Quinn model p0| , then the rate 
for ip — > 07 is proportional to x 2 and that for T — > cry is proportional to x~ 2 . The present bounds on these processes 
pof are only consistent with values of x ~ O(l)- Thus short-lived axion models to be viable must have variant quark 
couplings, |14[| such that both the ip — * 07 and T — * 07 rates are either proportional to x 2 or x~ 2 . 

These variant models are ruled out by a combination of other experiments, jljj The process n + — * e + e~e + i> e 
observed at SIN jl8| with a branching ratio of O(10 -9 ) bounds the parameter A3 in these models, since the decay 
chain n + — > ae + i/ e followed by a — > e + e _ would contribute to the signal. One finds |Q from these considerations 
that A 3 < 0.2. On the other hand, the I = to I = Ml decay of the 3.58 MeV excited state of 10 B to the ground 
state 19 1 provides a bound on the isoscalar mixing parameter A G A Q < 2. Although there are models in which 



5 This requires that, effectively, only the first generation of quarks (N g = 1) feel the C/(1)pq symmetry. 
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one can satisfy either one, or the other, of these two constraints, one cannot satisfy both because these parameters 
are themselves constrained by the relation 0] 

m a v AmT 

A 3 -A ~ — , — > 15 . (56) 

m* fit V ™d 



INVISIBLE AXION MODELS 

Although it was a sensible assumption to suppose that the J7(1)pq breaking scale / was the same as the weak 
scale v, this is not necessary. The dynamical adjustment of the strong CP angle — ► works for any scale /. If 
f v, the resulting axions are very light (m a ~ 1//), very weakly coupled (coupling ~ 1//) and very long lived 
(t(<z — > 27) ~ / 5 ) . Thus these axions are, apparently, invisible. 

Because / 3> v by assumption, in invisible axion models the C/(1)pq symmetry must be broken by an SU(2) x Z7(l) 
singlet VEV. Hence these models all introduce an 577(2) x f7(l) singlet complex scalar field cr. The invisible axion is 
then, essentially, the phase of a. Hence, concentrating again only on the axion degrees of freedom, one can write 

a = 4=e ia// . (57) 

Broadly speaking, one can classify invisible axion models into two types depending on whether or not they have 
direct couplings to leptons. The, so called, KSVZ axions |^0| are hadronic axions with only induced coupling to 
leptons. The, so called, DFSZ axions pi| ) on the other hand arise in models where axions naturally couple to leptons 
already at tree level. I describe below these two types of invisible axions in some more detail. 

KSVZ Axions 

In these models one assumes that the ordinary quarks and leptons are PQ singlets. The SU(2) x J7(l) singlet field 
cr, however, interacts with some new heavy quarks X which carry f7(l)pq charge |2C[] via the Yukawa interaction: 



-KSVZ 



-hX h <rX R - h*X R <r^X L . (58) 



The interactions of the KSVZ axions with ordinary quarks arises as a result of the chiral anomaly which induces a 
coupling 



r - - 

^anomaly « 



where qx is the electric charge of the heavy quarks X. Calculations analogous to the ones discussed earlier give the 
following axion parameters for the KSVZ model B 

\ 1 \ 1 (m d - m u ) 1 

A m = l; A 3 = ; A c = -- (60) 

^ (md + rn u ) Z 



and 



2 4md + m u 

K ail = 3q x - — . (61) 

3(m d + m u ) 



DFSZ Axions 

In this class of models plj , both the quarks and leptons carry ?7(1)pq and hence one again needs two Higgs fields 
$1 and $2- However, the quarks and leptons feel the effects of the axions only through the interactions that the field 
a has with $1 and $2 in the Higgs potential. This interaction, conventionally, occurs through the term 

£axion = «$f C$ 2 (a t ) 2 + h.c. , (62) 
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which serves to fix the PQ properties of a relative to $1 and $2. 

The contribution of the axion field a in $1 and $2 is readily isolated 0. Defining, as before, v 2 + v\ — v 2 , one has 



vi 



2v|a 



— exp I *— a 



$2 = 



V2 

V2_ 

71 



2v\a 
/ 



exp z 



(63) 



(64) 



where X\ = 2v 2 /v 2 ; X2 = 2v\jv 2 . These formulas correspond to what was done earlier (cf. Eq. (39)) with the 
replacements: x X\\ x _1 <-> X2. Using instead of / the scale 



/ = f/2N g 



(65) 



in the formula for the axion mass and mixing parameters (Eqs. (35)-(37)) gives the same expression for the axion 
mass in the DFSZ and KSVZ models. Adopting this convention, the axion parameters for the DFSZ case are 



A 
A 



1 



Aa 



Xi - X 2 _ (m d - m u ) 
(m d + m u ) 



2N n 



and 



2N n 



4m„ 



3 3(m d + m u ) 
Note that the DFSZ axion has a coupling to electrons given by 



C 



DFSZ 



Xi m e _ 
-i^ttz — —aej^e . 



2N g f 



(66) 



(67) 



(68) 



ASTROPHYSICAL AND COSMOLOGICAL CONSTRAINTS ON INVISIBLE AXIONS 

Invisible axion models are constrained substantially by astrophysical and cosmological considerations, which restrict 
the allowed range for /(/) — or equivalently the range for the axion mass 



m a ~ 6.3 



10 6 GeV 



/ 



eV . (69) 



The astrophysical bounds on m a arise, essentially, because axion emission removes energy from stars altering their 
evolution. Because the axion couplings to matter are inversely proportional to /, for large values of / (and therefore 
small axion masses) axions are less effective at cooling stars. Hence astrophysics provides a lower bound for / or, 
equivalently, an upper bound on m a . If axions are sufficiently light, their emission eventually is irrelevant for the 
evolution of stars. 

The bounds on invisible axions one derives from stellar evolution considerations are less restrictive for KSVZ axions. 
It turns out that the most effective cooling process is Compton axion production p2| je — > ea, which is proportiional 
to the coupling of axions to electrons. This process is absent for KSVZ axions, since gf^ /z = 0. However, one can 
obtain a bound on these axions since they can remove energy from stars via the Primakoff process 7Z, e — > Z,ea 
involving the axion to 2-photon coupling K a ~ n . 

Schematically, the energy lost by a star from the Compton process is given by 



1 

Pstar 



dn e dn-y\v\aE a , (70) 
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where p sta r is the stellar density and one integrates the interaction rate \v\a weighted by the axion energy over the 
density of the initial states. For m a <C T star ~ O(KeV), the only dependence on m a are through the dependence 
of the interaction rate on /: |u|er ~ f~ 2 ~ m 2 ,. Bounds are then obtained by requiring that Q < Qnucrthe rate of 
nuclear energy generation (Qmici ~ 10 2 ergs/g sec). 

From a detailed and careful study of how axion emission would affect stellar evolution, Raffclt |22] gives the following 
bounds on the two classes of invisible axions 

10~ 2 0.27 
(m ) D FSz < — r^eV ; (m a ) K svz < -T7 — eV . (71) 

Stronger bounds than these can be derived from the observation of neutrinos from SN 1987a. The bounds arise 
because if the axion luminosity is comparable to the neutrino luminosity (~ 10 53 ergs/sec) during the core collapse, 
then the neutrino signal would be altered. It turns out that the dominant process for axion production during the 
collapse is axion bremstrahlung off nucleons (N + N^N + N + a). As a result the SN 1987a bounds one obtains 
are quite similar for KSVZ and DFSZ axions. The results that Turner [p3j gives from his study of this issue are 

, . 1.7 x lO- 3 Tr 
(m a ) D FSz < — 72 eV ; 

^brcm 

(™ a ) K svz < 8.4 x 1CT 4 eV (72) 

where 

eL m = l-44+i(X 1 -X 2 - 1.55) 2 . (73) 

Cosmology, on the other hand, provides an upper bound for / (/) or, equivalently, a lower bound for the axion mass. 
This bound was derived by a number of authors p4| and its origin is easy to understand. When the Universe goes 
through the £/(1)pq phase transition, at temperatures of order T ~ /, the axion field acquires a vacuum expectation. 
At this stage the color anomaly is not effective, so the axion is a Nambu-Goldstone boson and (a p h ys .) ~ /• As the 
Universe cools to a temperature T* ~ Aqcd, the axion gets a mass of order m a ~ Aq CD // and the axion VEV is 
driven to zero dynamically (a p h ys .) — * 0, corresponding to 9 = 0. The relaxation of (a p h ys .) to this value is oscillatory 
and this coherent oscillation of the zero-momentum component of the axion field contributes to the Universe's energy 
density. The larger / is, the larger the axion contribution to the energy density of the Universe is. Asking that this 
contribution not exceed the Universe's closure density then gives an upper bound on /: 

A little more quantitatively Q, one can examine the equation of motion for (a p h ys .) in the expanding Universe in 
the approximation that only the axion mass term is relevant in the axion potential: 

d 2 (a p hy S .) R(t) d(a phys .) 2 

dp + 3 W)^^~ ra » W(aphys - ) - ' (74) 

where R(t) is the cosmic scale parameter and R(t)/R(t) = H(t) is the Hubble constant. If H(t) ^> m a (t) then there 
are no oscillations, while in the reverse limit the oscillations are sinusoidal. Oscillations start at T*, when 

A 2 

m a (T*) ~ H(T*) ~ — . (75) 

AlPlanck 

At the start of oscillations, the energy density in the axion field is of order p a (T*) ~ m 2 (T*)/ 2 . If one assumes that 
m a (t) is slowly varying then one can show that p4| p a (t) ~ m a (t)/R 3 (t). Thus the contribution of axion oscillations 
to the Universe energy density today is of order 



Pa = Pa{T*) 



m a (T*) 



R 3 (T* 
R?~ 



A y° T3 , (76) 



Planck 



where T ~ 3°K is the temperature of the Universe today. Requiring that p a be less than the closure density of the 
Universe provides the lower bound on m a . Amazingly, the above order of magnitude formula gives the same bounds 
on m a that the more careful calculations give Q 

m a > (10~ 5 - 1(T 6 ) eV . (77) 
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LOOKING FOR INVISIBLE AXIONS 



The cosmological and astrophysical constraints on invisible axions allow only a rather narrow window for the axion 
mass, roughly 

1CT 6 eV < m a < 10~ 3 eV . (78) 

If the axion has a mass near the lower limit, then axions play a very important role in the Universe and could be 
(part of) the Universe's cold dark matter. The cosmological bounds are subject to caveats, so it is useful to explore 
all masses below 10~ 3 eV. For instance, decay of axionic strings |25[| may increase the contribution of axions to 
the present Universe's energy density by about a factor of 10 3 -thereby pushing the axion mass lower bound to the 
astrophysical limit. On the other hand, Linde |26| has argued that in inflationary models, it is quite possible that 
(a p h ys (/)) ~ cf 1 with c <C 1. This would decrease the lower bound on the axion mass by a factor of c 2 . 

These uncertainties notwithstanding, it is very worthwhile trying to detect a possible signal of invisible axions. 
Indeed, experiments are presently under way to try to detect such axions on the assumption that they constitute the 
dark matter component of our galactic halo. Under this assumption one can estimate both the presumed axion energy 
density, p£ al ° - 5 x lCT 25 g cm 3 ~ 300 MeV cm , and their velocity, v a — 10 3 -the virial velocity in the galaxy. 



The basic idea for these experiments was suggested more than a decade ago by Sikivie |27|. It makes use of the axion 
to 2-photon coupling to resonantly convert in a laboratory magnetic field the halo axions into photons, which can 
then be detected in a cavity. 

The interaction of halo axions with a constant magnetic field B in a resonant cavity, as a result of Eq. (37), 
produces an electric field with frequency uj = m a . The generated electromagnetic energy can be detected in the 
cavity. When the cavity is tuned to the axion frequency, one should observe a narrow line on top of the noise 
spectrum. On resonance, one can write for the axion to photon conversion power the expression ps| 



„halo 
Pa 



[VBl] 



— K 

f -"-077 



CoverQoff • (79) 



Here the first factor is just the number of axions per unit volume; the second details the magnetic energy; the third 
is the axion coupling strength (g 2 77 ); the fourth is an overlap factor of O(l) which depends in detail on which mode 
is excited in the cavity. Finally the last factor is the effective Q for the experiment, which is the minimum of the 
cavity's Q and that produced by the spread in the axion frequency due to its velocity. Typically both numbers are of 
O(10 6 ). 

Halo axions produce microwave photons (1 GHz = 4 x 10~ 6 eV). Two pilot experiments, done at Brookhaven 
p9| and the University of Florida had limited "magnetic energy" -typically B 2 V ~ 0.5 (Tesla) 2 m 3 -and relatively 
noisy amplifiers. These experiments set limits for g 2 77 about a factor of 100-1000 away from theoretical expectations. 
There are now, however, two second generation experiments which should be sensitive at, or near, the theoretical 
limits. 

The first of these experiments, which has just started taking data, is being carried out at Livermore National Lab 
pl| . It has a very large magnetic volume (B a V 2 ~ 12 (Tesla) 2 m 3 ) and state of the art amplifiers to reduce the noise 
levefl The second experiment, to be carried out in Kyoto University [p2j , uses a rather moderate magnetic volume 
(B Q V 2 ~ 0.2 (Tesla) 2 m 3 ). However, it employs an extremely clever technique for counting the produced photons 
using Rydberg atoms, which is extremely sensitive. I show in the Figure the expectations of these experiments, along 
with the limits obtained by the pilot experiments. It is clear that these experiments have the capability to answer 
the question of the existence of axions in the mass range that they can probe. Results are expected in the next few 
years. So, stay tuned! 
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Even so the power produced by halo axions is tiny. At 1 GHz one expects Paxion ~ 4 x 10 Watts! 
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Experiments done by: 




Axion Mass [micro-eV] 



Exclusion regions in the axion mass versus coupling constant plane obtained by the pilot cavity experiments. Also shown are the 
coupling constants expected in a range (DFSZ and hadronic) of axion models. The area extending into the KSVZ (hadronic 
axion) region is the expected sensitivity of the U.S. experiment. The narrow region extending into the SFSZ region is the expected 
sensitivity of the Japanese experiment, (See Ref. [2SJ.) 



